Impurity-induced polar Kerr effect in a chiral p-wave superconductor 
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We discuss the polar Kerr effect (PKE) in a chiral p-wave (p x + ip^-wave) superconductor. It 
is found that the off-diagonal component of a current-current correlation function is induced by 
impurity scattering in the chiral p-wave condensate, and a nonzero Hall conductivity is obtained 
using the Kubo formula. We estimate the Kerr rotation angle by using this impurity-induced Hall 
conductivity and compare it with experimental results [Jing Xia et al., Phys. Rev. Lett. 97, 167002 
(2006)]. 
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Recently, the quasi-two-dimensional (quasi-2D) super- 
conductor Sr2RuC>4 with T c = 1.5 K has attracted 
considerable attention and it has been investigated 
extensively^ It is plausible that the order parameter in 
Si'2Ru04 has the spin-triplet p x ± ip y -~w&ve symmetry. 
One of the fascinating properties of this state is the 
spontaneous breaking of parity in a 2D sense (p x — > p x , 
p y — > —p y ) and time-reversal symmetry due to the pres- 
ence of nonzero chirality characterized by l z = ±1, where 
l z is the z-component of the relative orbital angular mo- 
mentum of the Cooper pair. 

The polar Kerr effect (PKE), in which the direction of 
polarization of reflected linearly polarized light is rotated, 
has been known as an effective tool for understanding 
ferromagnetism.— Because of the analogy between ferro- 
magnetic order, for instance, with s z = 1 and chiral pair 
condensation with l z = 1, it is naively expected that the 
PKE is induced in the chiral p-wave state at zero field. 
In fact, the PKE has been observed in the superconduct- 
ing state of Sr2RuC>4<^ Up-to-date theoretical reports on 
the PKE in the chiral p-wave state are given in Ref4^, 
in which interesting mechanisms have been proposed by 
the field theoretical approach; however, obtained results 
of the Kerr rotation angle are considerably smaller than 
the experimental results. 

Therefore, it is crucial to elucidate the fundamental na- 
ture of the PKE in the chiral p-wave superconductor. We 
will show that a Kerr rotation angle comparable to that 
obtained experimentally is obtained by taking into ac- 
count nonmagnetic and short-ranged impurity scattering 
of quasiparticles in a chiral p-wave condensate. It is also 
found that this impurity-induced PKE is suppressed or 
zero for any superconducting state other than the chiral 
p-wave stated This result is contrary to the naive analogy 
with a ferromagnet, since the effect is not proportional to 
chirality analogous to magnetization, but suppressed in 
higher-chirality states with l z = ±2, ±3, ■ • •. The natural 
unit h = c = ks = 1 is used throughout this study. 

We review the phenomenology of the PKE in time- 
reversal-symmetry-breaking superconducting systems, 
which is an extension of the argument for itinerant fer- 
romagnetic systems.— & We also refer to the discussion of 
anyon superconductivity^ Suppose that z > is empty 
and z < is filled by the superconductor and inci- 



dent light is linearly polarized and propagating along the 
^-direction perpendicular to the superconducting plane 
with a wavevector q = — zq z . The Maxwell equa- 
tions inside the material are z\7 z x E = —dH/dt and 
zV 2 xB=j + d~E/dt. It is obvious that for manifestly 
gauge invariant systems, current ji = J2j= x y a ij^j — 

x,y, o~ij is the con- 
ductivity tensor and A4 and <j> denote the vector and 
scalar potentials, respectively. In superconductors, 



J2j yAj + d<p/dxjj , where i 



and in general, erjj-' ^ af- 1 because of the spontaneous 
breaking of the gauge symmetry.—^ The scalar potential 
becomes redundant in this problem and we select a <p = 
gauge. In this gauge, the Maxwell equation inside the 
superconductor is 



(1) 



I (j 2 — q\ + ibso 



Sij + zce^M} A M (w) = 0; 



here, we suppose that the long- wavelength limit q << 
£ , where £ is the coherence length of the supercon- 
ducting order parameter and we omit the q dependence 
of the conductivity tensor. It is clear that there are two 

propagating modes qf — \J uo 2 + iujg x v x (ui) ± Loo x V y (ui) . 

In Refi^, the low- frequency limit w « 2|A| has been 

considered, and only the static values cjj ( w = 0) 
have been considered, however, this value is not suit- 
able for the experimental situation in the Ruthenate 
u = 0.8 eV » 2|A| ~ 1CT 4 eV Following Ref4, we 
solve the Maxwell equation with an appropriate bound- 
ary condition at z = and obtain the Kerr rotation angle 



e K = -Im 



- q z ) 



qi q? 



(2) 



which would be applicable in a wide frequency region. 
The factor (q+ — q~) indicates that a x v y \uj) is crucial to 
the PKE. 

Let us calculate a^\uj) with lu » 2|A| in the chi- 
ral p-wave state. We use a quasiparticle Hamiltonian 
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FIG. 1: Leading diagrams of eij7rij(iu n ) . The lines 
with the arrows —>—, — < — >—, and — > — < — 
denote quasiparticle Green functions G p (ie m ), 3~p(ie m ), and 
J-p(ie m ), respectively, in G p (ie m ). The lower lines have fre- 
quency e m +„ = e m + uj n , while the upper lines have e m . The 
dashed lines denote the impurity coupling v\. The scattering 
at the black dots yields the factor Tr(r3g p ■ rgk • t) in Eq. 
©■ 

with a cylindrical Fermi surface that models the dom- 
inant 7-band in Sr 2 Ru0 4 i The elect romagnetic inter- 
action and impurity scattering are also taken into ac- 
count. In the Nambu representation ^ p = (c p j, c*_ p ^) T , 
the Hamiltonian H = Ho+H em +Hi. The first part Hq = 

I (€f*p§p ' T *P' where gp = (&Ap,-ImA p ,e p = 

2 

s^— — cp), t — (ti, T2, t 3 ) is the Pauli matrix in the 
Nambu space, and A p = \A\(p x +ip y ) is the momentum- 
dependent part of a chiral p-wave gap function with p = 

p/|p|. The second part H em = J (j<f> + j^) -A q is 
the electromagnetic interaction with the Fourier forms of 
the paramagnetic current ffi = -e / (5^* P ^7* P + q 

=~ e I (2?o :j ^p^p.p+q^p+q an< ^ the diamagnetic current 
jq "i/lCF^^p+qAq. The last part H { = 
v ' 1 1 (2t!\ 3 <°iqi°-q i s ^ ne nonmagnetic impurity scattering 



with the Fourier components of the quasiparticle den- 
sity pq = J ^^.^prs^p+q and impurity density p iq = 
^ giq-Rj impurity site). It is assumed that the 

scattering potential is short ranged and the s-wave chan- 
nel is dominant. 

Let 7Ty(w) denote the Fourier component of the two- 
point retarded correlation of the paramagnetic current 
jqXo- The Hall conducvitity is obained by the Kubo 

formula as^ <7xy(w) = -K^ y {uj)/ioj = e y ,7r y( w )/ 2 ^) 
where eij is the totally antisymmetric tensor in the 2D 
plane. To obtain Ylij e ij 7r ij( u ')'> we calculate the Mat- 
subara form . eijirij(iui n ), where u> n = 2nirT is the 
bosonic Matsubara frequency at temperature T, and use 
the relation ^"^M = J2ij ^ij^ij^n —> w + £7), 
where 7 denotes dissipation*^ The Matsubara Green 
function for quasiparticles is given by 

r ( U \ _ 1 = ( Gp( ie rn) Fp{iem) \ 

p[ m) ~ ie m +E P -T-{f;(-iem) -0 P (-iem) J ' 

where e m = (2m + l)irT. 

In literature ; 11 ' 12 ' 13 ' 14 ' 15 it has been pointed out that 
there is no contribution of the one-loop diagram to the 
Matsubara form eijiTij{iu) n ), i.e., the zeroth-order 
term of the impurity scattering since 7i, P , P = Pi/m e . 
Then, vertex corrections must be taken into account to 
obtain nonzero contributions. The leading contribution 
is shown in Fig. 1, which is in the first order of the 
impurity concentration n; and in the third order of the 
impurity strength v x . These diagrams are similar to the 
skew scattering diagrams in the extrinsic anomalous Hall 
effect iii We obtain 



E 



eij7Tij(iu} n ) 
T 



iitfe 2 7/3 ^ {Tr(Gp(ie„j7 l! p i pGp(ie m+ „)T 3 Gk(«em+r07j,k,kGk(iem)T3Gk 1 («em)T-3) 



mpkki ij 

)r 3 G k ( l )7j,k,kGk(ie m )r 3 )} + 0(n\vf) 



T 



V 3 ^ 



n\vfe 2 v 2 F i> x kTr (r 3 g p • xg k • r) u n (e m + e m+n ) 2 (e m e m+n - E^J 



(3) 



where E p = |g p | = y'e 2 + |A p | 2 and e m+n = e m + ui n . 

We use the cylindrical coordinate for representing mo- 
mentum integration. As observed from the leading term 
of Eq. ([3]) , an azimuthal dependence p x k — sin 9 p k (9 p k ■ 
angle between p and k) arises from the contraction of the 



vertex part £V e^.p^.k.k, and 
Tr (r 3gp • Tg k • t) = -2i (ReA p ImA k - (p <-► k)) (4) 

arises from the impurity scattering at the black dots 
in Fig. [T] In the chiral p-wave (p x + ip y -wetve) state, 
Tr (r 3 g p • rgk • t) = 2i|A| 2 sin^pfe. Then, the leading 
term of Eq. ([3]) survives after the azimuthal integration. 
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It is obvious that this term yields the second-order contri- 
bution of the gap amplitude |A(T)| 2 - |A(0) | 2 (1 -T/T c ) 
in the Ginzburg-Landau (GL) regime. By using GL ex- 
pansion, we estimate the large uj n part of the term to per- 
form momentum integral first and Matsubara sum later. 
By using analytic continuation and the Kubo formula, 
we obtain 

<V) = ilcs (i-pir (—) 3/2 7 f, , 3 > ( 5 ) 

V T c J£q\itt q J (uj + i/tq) 6 

where a x °J = e 2 /2ird (d: layer distance), k = (njd) -1 / 2 , 
t _1 = n i vfN(Q)/2 (N(0) = m/2-Kd: density of states at 
the Fermi surface per spin), £o = Vf/ttT c is the super- 
conducting coherence length at T — 0, £f is the Fermi 
energy, and j BCS = |A(0)|/T c ~ 1.8. 

In general, when we consider chiral states with l z = 
±1, ±2, • ■ •, i.e., (p x + ip y ) lz -wa,ve states, Eq. ((4]) be- 
comes 2i\ A| 2 sml z 8 p k, and then, the leading term of Eq. 
([3]) vanishes after the azimuthal integration, except for 
l z = ±1, i.e., the chiral p- wave state. If we consider chiral 
states with a horizontal line node such as (p x + ip y ) lz p z , 
Eq. (UJ) becomes 2i\ A\ 2 p z k z sin l z 9 p k and the leading 
term also vanishes for any l z after integrating out the 
z-component of the momentums. For nonchiral (l z = 0) 
and time-reversal-breaking states, i.e., d + is-wave pair- 
ing, Eq. (U]) is proportional to cos2# p — cos2(?fc and 
the leading term will vanish after the azimuthal integra- 
tion. For nonchiral and time-reversal-symmetric states 
(s-wave, d x 2_ v 2-wave, ... ), a x v y \uj) should be zero be- 
cause of the symmetrical reason. In fact, it is easy to 
confirm that Eq. ([4]) becomes zero since the gap func- 
tion can be made real for these states by using the U(l) 
phase transformation of fermion fields. To sum up, in this 
impurity scattering mechanism, <j x V y (uj) is suppressed or 
zero for any state other than the chiral p-wave state (see 
Footnote^). 

The diagonal component a x v J (oj) is effectively approx- 
imated by the Drude form a x v J (uj) ~ w 2 To(l — iujtq)~ 1 
in the high-frequency limit uj >> 2|Ap^, where uj p = 
\f n e e 2 / m e is the plasma frequency and n e is the elec- 
tron number density. In fact, this behavior of a x v J(uj) 
has been verified experimentally in the superconducting 
state of Sr 2 Ru0 4 ^ 

Let us estimate Ok from Eq. @. We use suitable 
parameters for the experiment, i.e., d — 6.8 A, £o = 
660 A, uj = 0.8 eV, uj p = 1.3 eV, e F = 0.14 eV, and 
Tq 1 = 6.6 x 10~ 5 eV £ The impurity-mean distance k 
is rather ambiguous, and we consider a variation 1000 — 
5000 AM The estimated value of K from Eq. fjjj) at 
T/T c = 0.8 is 6 nanorad for k = 1000 A and 30 nanorad 
for li = 5000 A. The latter case agrees well with the 
measurement result of approximately 60 nanorad at the 
same temperature^ 

The frequency dependence of 0k in Eq. (J2J) at T/T c = 
0.8 is plotted in Fig. [5] We find that 0k behaves as uj^ 1 
for uj < uj p , as indicated by Eqs. §5§ and ©■ At u> = uj p , 
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FIG. 2: uj dependence of 9 K in Eq. © at T/T c = 0.8, ob- 
tained using parameters d = 6.8 A, h = 10000 A, £o = 660 A, 
cop = 1.3 eV, e F = 0.14 eV, and r _1 = 6.6 x 10 -5 eVA It 
is assumed that the cut off energy of the pairing interaction 
ujd > Wojp. 

0k shows a sharp peak and behaves as u>~ 3 for uj > lu p . 

The temperature dependence of 0k can be explained 
as follows. Eq. |(3J) is effectively approximated by 

K ~^ife(a#M), (6) 

UJ p 

when uj 2 — uj 2 » \wa^ v \uf)\ and ujtq >> 1 (these condi- 
tions are satisfied using suitable parameters in the exper- 
iment). The approximated form Eq. ^ indicates that 
0k depends linearly on the temperature in the GL region 
(see also Eq. (|SJ)). This result is not in agreement with 
the experimental report;^ however, the error bars of the 
data are large and it appears to be difficult to discuss the 
temperature dependence. 

The obtained result might be overestimated since it is 
assumed that every scatterer has the same fixed poten- 
tial strength rather than a distribution. We also suppose 
that the energy scale of the pairing interaction ujd is 
well above the laser frequency uj ~ 0.8 eV. Although 
the nature of the pairing interaction is not elucidated, 
ujd is likely to be less than the used laser frequency. It 
is pointed out that the Hall conductivity would be sup- 
pressed for uj > u>d£- However, we believe that the pre- 
sented mechanism plays an important role in the PKE in 
Sr 2 Ru0 4 £ 

We comment on previously proposed theories of the 
PKE in the chiral p-wave state. In Reb^i, a x S y in 

Eq. {!]) is obtained using the relation a X y(u>) = 
|e i j^-7r ( ^ (a;,q)| q= o, where 7r^(o;, q) is the correlation 
between charge density and current density. The calcula- 
tion has been performed in the case without any impurity, 

(s) 

and nonzero values of a xy have been obtained; however, 
as shown in Eq. @, a^y (uj) and not a x S y (uj) is responsi- 
ble for the PKE. In Refill, it is implicitly assumed that 
cr x S y (uj) = <r x V y (uj) in the high-frequency limit (see the dis- 
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cussion below Eq. (17) in Ref*2i. The same assumption 
has been made in Refi22); however, this would not be true 

since <Txy (to) is zero without impurity scattering for an 
arbitrary frequency . 11 ' 12 ! 13 ' 14 This problem has also been 
pointed out in recent arguments^ In Refs^£, aiy (uj, q) 
has been discussed. A remarkable finding is that the fi- 
nite q effect is responsible for the PKE (a similar argu- 
ment has been provided in Ref. 23 ); however, the obtained 
Kerr angle is about 9 orders of magnitude smaller than 
that obtained experimentally^ 5 - 

Let us turn to discuss the Chern-Simons term 
p, v = 0,1,2) which is generally in- 
duced in the low-energy and long-wavelength effective 
Lagrangian for gauge fields obtained by integrating out 
2D electrons in a system with parity and time-reversal- 
symmetry breaking^ The origin of the term is closely re- 
lated to the parity anomaly . 25 ' 26 ' 27 ' 28 It has been pointed 
out that for the chiral p-wave superconductor, a nonzero 
value of <Jxy (u) = q = 0) is obtained without im- 
purity scattering. At T = 0, a [ x s J{0) = fijll^^! 4 . 
the finite-temperature effect was investigated in the GL 
schem o 14 ' 29 (see also Ref.—). This indicates that a part 

of the Chern-Simons term — (A Q diAj + AidjA ) 
is induced in the effective Lagrangiam 11 ' 12 We can con- 
clude that this part has an "intrinsic" origin since it is 
induced without impurity scattering, while a x v y ; (0) = 0, 
and the other part of the Chern-Simons term eijAidoAj 

is not induced i 11 ' 12 ' 13 ' 14 As we have emphasized, o^y (w) 



becomes nonzero when we take into account vertex cor- 
rections from impurities, shown in Fig. [T] From the 
leading term of Eq. ((3J, we can observe that in the static 
limit, the second-order contributions of the gap ampli- 
tude vanish occasionally, and the forth-order contribu- 
tions give 



r(°) 



127C(7) k ( e F 



128tt 4 £o \nT T* 



T 



~4 

iBCS 



near T< 

<#>) 



This indicates that the other part 
eijAidoAj has an "extirisic" origin. 



In summary, we have discussed the PKE in a chiral 
p-wave superconductor. We have estimated the contri- 
bution of impurity scattering in the chiral p-wave con- 
densate. The skew-scattering-type diagrams (see Fig. [T]) 
show the leading contribution. In this impurity-induced 
mechanism, the PKE would be suppressed or zero for 
any state other than the chiral p-wave state (the possi- 
bility of nonunitary pairing is not taken into account (see 
Footnote 6 -). 
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